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Abstract 

We study quark flavour violation (QFV) in the squark sector of the Minimal 
Supersymmetric Standard Model (MSSM). We assume mixing between the second 
and the third squark generation i.e. cl,r — tL,R mixing. We focus on QFV effects 
in bosonic squark decays, in particular in the decay into the lightest Higgs boson, 
v,2 uih°, where ui^ are the lightest up-type squarks. We show that the branch- 
ing ratio of this QFV decay can be quite large (up to 60 %) due to large QFV 
trilinear couplings, and large c~l,r — ti,R and £l — £r mixing, despite the strong con- 
straints on QFV from B meson data. This can result in remarkable QFV signatures 
with significant rates at LHC, such as pp — > ggX — > u\^iu2cX — > ui^iui^cX — > 
cxiicx±h°cX(= cicch°]^TX). The QFV bosonic squark decays can have an influence 
on the squark and gluino searches at LHC and can play a role in the determination 
of the MSSM parameters, in particular of the QFV trilinear couplings. 



1 Introduction 



In most searches for supersymmetric (SUSY) particles at the LHC, the analyses have been 
performed within simplified SUSY models. However, SUSY extensions of the Standard 
Model (SM) can have a richer structure. In the Minimal Supersymmetric Standard Model 
(MSSM), mixing between the different squark generations is in principle possible. This 
could lead to quark flavour violating (QFV) effects, in addition to those induced by the 
Cabibbo-Kobayashi-Maskawa (CKM) matrix [1, 2, 3]. Mixing between the 1 st and the 2 nd 
squark generation is strongly supressed by K physics data [4]. Therefore, in this paper we 
assume mixing between the 2 nd and the 3 rd squark generation, respecting the constraints 
from B physics (see Appendix B). Although these constraints are also quite severe, they 
allow nevertheless substantial QFV effects. 

In the MSSM, the mixing of the 2 nd and the 3 rd squark generation was theoretically 
studied for squark and gluino production and their decays at the LHC in the context of 
Minimal Flavour Violation (MFV) [5, 6, 7] as well as for general flavour mixing [8, 9, 
10, 11, 12, 13, 14, 15]. As shown in these papers the effects of QFV can be large. For 
example, in the case of mixing between top and charm squarks, we can expect a large 
branching ratio (up to 40 %) of the QFV decay of the gluino, g — > ci(ct)x\ [13]. This is 
due to the fact that the lightest up-squark mass eigenstates Ui >2 are mainly mixtures of 
tfi and cr. Hence, iti and u 2 can both decay into CX\ and txi- 

In addition to the fermionic decays of squarks there are bosonic decays, qi — > cjj + 
Z ,h°,H°,A and q^ — > q'- + W ± ,H ± , if kinematically allowed. In the quark flavour 

conserving (QFC) case, the most interesting decays are t 2 — > t± + Z°, h°, H°, A ; b 2 — > 
h + Z°, h°, H°, A ; b 2 -> h + W~ , H~. They were studied in [16, 17]. The QFV bosonic 
decays were recently considered in [11]. There the characteristic differences to the MFV 
case were worked out. A non-minimal flavour structure in the squark sector can change 
the entire squark decay pattern quite drastically, because many more transitions are 
possible. 

In the present paper we are particularly interested in the bosonic QFV squark decays 
into the lightest Higgs boson, qi — > q~j + h° in the MSSM. These decays offer the best 
possibility of determining the trilinear couplings q^ — qj — h° entering the soft-SUSY- 
breaking Lagrangian. Another possibility would be to study the 3-body production pp — >■ 
q~iq~jh as discussed for example in [18] for the QFC case. As the Higgs boson couples 
dominantly to the qh~qR combination, one gets information from the decays q~i — > q~j + h° 
on the flavour structure of the left-right (LR) terms in the squark mass matrix. We shall 
study the mixing between the 2 nd and the 3 rd generation of up-type squarks i.e. i L ^ R — c LtR 
mixing. There are strong constraints on this mixing from B physics, Higgs boson searches 
and SUSY particle searches (see Appendix B). In our analysis we take into account all 
these constraints. 

The paper is organized as follows: In Section 2 we shortly give the definitions of the 
QFV squark mixing parameters. In Section 3 we discuss the QFV bosonic decays of up- 
type squarks in detail in a definite scenario. We also consider two further scenarios, one 
GUT inspired and another one where the bosonic decays of u 2 dominate over the fermionic 
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decays. Section 4 contains a discussion of various QFV signatures. In Section 5 we give 
a summary. In the Appendices we show explicitly the part of the interaction Lagrangian 
we use in our analysis and summarize the experimental and theoretical constraints on the 
MSSM parameters, especially those on the QFV parameters mainly from B physics. 



2 Squark mixing with flavour violation 

In the MSSM the most general form of the squark mass matrices in the super-CKM basis 
of g 07 = (?iL,g2L,g3L,giR,g2R,g3R), 7 = 1,-6, with (91,92,93) = (u,c,t),(d,s,b) is [19] 

" \ M\ RL M\ RR )• U 

for q = u,d, where the 3x3 matrices read 

M \ll = M i + D J,ll1 + rhl M\ LL = V CKM M 2 Q V^ M + D^ L 1 + m 2 , 

M Irr = M l + DdfiR 1 + m l M lRR = M u + Du,rr1 + ml (2) 

Here Mq^u,d are the hermitian soft-SUSY-breaking mass matrices of the squarks and 
m Ut d are the diagonal mass matrices of the up-type and down-type quarks. -D<j,ll = 
cos2/3m|(T| — e^sm 2 ^) an d ^ = e q sin 2 6w cos 2(3m 2 z , where T3 and e q are the 
isospin and electric charge of the quarks (squarks), respectively, and 6 W is the weak 
mixing angle. The left-left blocks of up-type and down-type squarks are related by the 
CKM matrix Vckm due to the SU(2) L symmetry. Note that VckmMqVq KM ~ Mq as 
Vckm — 1- The off-diagonal blocks of eq. (1) read 

M Irl = M Ilr = ^-^tan/3, 

Ml RL = M 2 l LR = ^=TZ-SrhuC0t(3, (3) 

Here Ty D are the transposes of the soft-SUSY-breaking trilinear coupling matrices of the 
up-type and down-type squarks Tu,d defined as Ci nt D —{Tuapu^ULaH^+Toap^pdLaH^), 
pi is the higgsino mass parameter, and tan/3 = i^/ui, where v\ : 2 = v^2(if° 2 ) are the 
vacuum expectation values of the neutral Higgs fields. The squark mass matrices are 
diagonalized by the 6x6 unitary matrices B?, q — u,d, such that 

B?M\{R«y = diag(m|, . . . , m|) with m- qx < ■ ■ ■ < m 9 ~ 6 . (4) 
The physical mass eigenstates q^i = 1, 6 are given by qi = Rj a qo a . 
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We define the QFV parameters in the up-type squark sector S^, ^ RR and 5% RL 
{a ^ ft) as follows [20]: 



*5 = K«8/\M*M*, (5) 



6 U J R = M^/^M* Uaa M^, (6) 

5 uRL ee (v2/V2)T ufja /^Mi aa M^ . (7) 

Here a, (3 = 1,2,3 (a ^ (3) denote the quark flavours. The QFV parameters relevant for 
this study are S^ L , S^ R = (6%? L )*, S^ R and 5^, which are the c R - t L , c L - t R , c R - t R 
and cl — tL mixing parameters, respectively. We also use the QFC parameter 5% RL which 
is defined by eq. (7) with a = (3 = 3 and is the II — t R mixing parameter. We assume all 
QFV parameters to be real. 



3 QFV bosonic decays of up-type squarks 

If kinematically allowed, the following QFV bosonic decays of up- type squarks are possible 
iH -> Uj + h°, H°, A (8) 

Ui ->■ dj + H + (9) 

in ->■ Uj + Z° (10) 

in ->■ dj + W + (11) 

with i,j = 1,...,6 specifying the squark mass eigenstates which are mixtures of the squark 
flavour eigenstates (see Section 2). Here h°(H°) is the lighter (heavier) CP-even neutral 
Higgs boson, A is the CP-odd neutral Higgs boson and H + is the charged Higgs boson. 
Of course, there are also QFC bosonic squark decays. In this article we study mainly 
""1,2,3 decays in scenarios where their decays into charged bosons of eqs. (9) and (11) 
and those into the heavier Higgs bosons H° and A are kinematically forbidden. The 
couplings between Ui — Uj/dj and the bosons in eqs. (8) - (11), taking into account QFV, 
are given in [11]. For completeness, the couplings to the lightest Higgs boson, h°, are listed 
in Appendix A. Note that the QFV parts are proportional to the soft-SUSY-breaking 
trilinear coupling parameter Tu- In the following discussion of the decays we will adopt the 
QFV parameters 5^, 5^ R ,S^ L ,S^ R (= 5^ L *) as defined in Section 2. The parameters 
S%3 R , S% RL are proportional to Tjj 23 and Tjj 32 , respectively. 

In the calculation of the branching ratios of the decays (8) - (11) we have to take into 
account both QFV and QFC fermionic squark decays [10, 13] 

Ui->u a + xl (12) 
Ui-^d a + xt (13) 
Ui^u a + g (14) 
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Table 1: Weak scale basic MSSM parameters at Q — 1 TeV for Scenario A, except for 
m A o which is the pole mass (i.e. the physical mass) of A . All of Tu aa and To aa are 
zero, except for T U33 = —2160 GeV (i.e. S^ L = —0.34). All other squark parameters not 
shown here are also zero. 



Mi 


M 2 


M 3 


A* 


tan j3 


rrijifi 


400 GeV 


800 GeV 


1000 GeV 


2640 GeV 


20 


1500 GeV 





a = 1 


a = 2 


a = 3 




(2400) 2 GeV 2 


(2360) 2 GeV 2 


(1450) 2 GeV 2 


*4aa 


(2380) 2 GeV 2 


(780) 2 GeV 2 


(750) 2 GeV 2 


M Daa 


(2380) 2 GeV 2 


(2340) 2 GeV 2 


(2300) 2 GeV 2 



XuLL 
°23 


zuRR 
°23 


zuRL 
°23 


xuLK 
°23 


0.024 


0.3 


-0.07 






Table 2: Physical masses in GeV of the particles in Scenario A (see Table 1). 





m x% 




A-4; 


Xi 




397 


824 


2623 


2625 


825 


2625 


















m h o 


m H o 




m H + 






126.0 


1496 


1500 


1510 





m- g 




m« 2 




m« 4 


mu 5 




1141 


605 


861 


1477 


2387 


2401 


2427 



where a — 1, 2, 3 is the flavour index, /c = 1, 4 are the neutralinos and xt ■> ^ = 1 ; 2 
are the charginos. As U\^. are mainly mixures of cr, and in the scenarios under 
consideration, both decays U\^. tXi an( i ^1,2 ~~ ^ C X? are possible. 

In the following we will discuss the QFV bosonic squark decays in three different 
scenarios. The first one which we call Scenario A and which we will study in detail 
is shown in Table 1. The corresponding basic SUSY parameters are given at the scale 
Q — 1 TeV according to the SPA convention [21], except for m^o being the pole mass of 
A . For this scenario all experimental and theoretical constraints given in Appendix B 
are satisfied, with the following values of the low energy observables: AMb s = 17.2 
ps" 1 , B(6 ->■ S7) = 2.9 • 10" 4 , B(B S ->■ /i+/i") = 3.3 • 10" 9 . All numerical calculations 
in this study, except for the cross sections, are performed with the public code SPheno 
v3.1 [22, 23]. We also use the package SSP [24]. Further on, we choose relatively large 
values of M 1; M 2 and /i in order to avoid the dominance of the fermionic squark decays. 
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Figure 1: Dependence of the masses of Ui (a) and u<i (b) on 5%% and in Scenario 
A, with the parameters given in Table 1. "X" in both plots corresponds to the reference 
point. 
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Figure 2: The mass of the lightest Higgs boson, m h o, as a function of T;y 33 and /i (a) and 
as a function of and 5^ (b) in Scenario A, with the parameters as given in Table 1. 
"X" corresponds to the reference point in both plots. 



M 3 is chosen such that the gluino mass is within the reach of LHC (14 TeV). Note that in 
the present scenario M 1; M 2 and M 3 do not fulfill gaugino mass unification at the GUT 
scale. The GUT-inspired case with gaugino mass unification will be studied in a second 
scenario. The physical masses of squarks, gluino, charginos, neutralinos and Higgs bosons 
are shown in Table 2. We obtain m^o = 126 GeV which is in the range of the possible 
Higgs signal at LHC [25] (see Appendix B). Moreover, in this scenario we are in the 



6 




(a) (b) 



Figure 3: The coupling \c{u 2 U\h°)\ 2 as a function of 5 2 RL and S 23 R (a) and as a function 
of S 2 RL and 5 33 L (b) in Scenario A of Table 1. The reference point is indicated by "X" 
in both plots. 



decoupling limit with m^o = 1500 GeV ^> m^o, and hence the lightest Higgs-boson h° is 
SM like. 

In our previous study on QFV fermionic decays of gluinos and squarks it turned out 
that QFV effects mainly depend on 5 2 RR and S 23 R , whereas the influence of the other 
QFV parameters was much weaker [10, 13, 14]. In the present paper we concentrate on 
the dependence of the QFV effects on S 2 RL and S 23 R , which enter the squark-squark-Higgs 
couplings. 

In Fig. 1(a) and 1(b) we show the mass contours of U\ and u 2 in the 5 2 RL — S 23 R 
plane. In all contour plots in this article the white regions satisfy all the experimental 
and theoretical constraints listed in Appendix B. In this scenario with large 5% RR , 8 3 R 
and S 2 RL , ui is mainly clIr — cr^—Il) mixture and u 2 is mainly a cr — £l) mixture. 
One can see a somewhat stronger dependence on 5 2 RL due to the sizable mass-splitting 
induced by the CR — ti mixing, which is a consequence of the chosen hierarchy within Mq 
and Ml. 

In Fig. 2(a) the lightest Higgs mass m^o is shown as a function of T\j 33 and \i. In 
order to obtain m^o ~ 126 GeV a large \Tu 33 \ entering the stop loop corrections to m" 
and a rather large \i entering the sbottom loop corrections through the term A^ — fi tan /3 
[26, 27, 28] are required. We have fi = 2640 GeV and T U33 = -2160 GeV in Scenario A. In 
Fig. 2(b) we show contours of m^o in the 8 23 R - 5 2 RL plane. Within the QFV parameter 
range shown the m^o varies by about 3 GeV due to the c admixture in the stop loops. 

Next we will study the decay u 2 ~ > uih° in more detail. In Fig. 3(a) we show the 
squared coupling \c(u 2 itih )\ 2 (see eq. (17) in Appendix A) as a function of 5 2 RL (~ T U32 ) 
and S 23 R (^ T U23 ). Dominant terms in c{u 2 u\hP) are those proportional to T U33 (= i R — 
ti — H 2 coupling) and Tu 32 (= c R — ti — H 2 coupling) (see Ci nt in Section 2) since Ui t2 
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are mainly mixtures of c R , t R and t^, and h° ~Re(i?2). Note here that in our scenario 
there is a large c R -t R mixing, S 2 RR = 0.3, a large t L -t R mixing, S 3 RL = —0.34. The term 
proportional to TJ723 (= t R ~ Cl — H 2 coupling) is rather small since the Cl component 
in uip is small mainly due to very large rric L (~ 2.4 TeV) (see Table 1). The decrease 
of the |c({t2Ui/i°)| 2 with increase of 8 2 RL is due to the fact that the contributions of 
these two couplings, T U33 and T U32 , to c(u 2 Uih°) interfere with each other constructively 
(destructively) in the case T U33 < and T U32 < (T U33 < and T U32 > 0). Again the 
dependence on 5 2 RL is much stronger than that on 5 23 R because the contribution of the 
coupling T(723 to c(u 2 U\h°) is rather small due to the small Cl component in u\ )2 in our 
scenario. Fig. 3(b) shows \c{u 2 u\h°)\ 2 as a function of the QFV parameter 5 23 L and the 
QFC parameter S 33 R with the latter being the t L — t R mixing parameter. Also here a 
strong variation of | c{u 2 U\ h°)\ 2 can be seen. It increases with increase of \5 33 L \ since the 
dominant contribution proportional to Tu 33 in c(u 2 uih°) increases because of increasing 
of both T[/33 and the ti component in ui j2 . The decrease of the |c(-U2Wi/i°)| 2 with increase 
of 8 2 RL is due to the same reason as in Fig. 3(a). 

In Fig. 4(a) we show the branching ratio of the decay u 2 — > Uih° as a function of the 
QFV parameters 8 23 L and S 23 R . This branching ratio can exceed 50%. Analogous to 
Fig. 3(a), the dependence on 5 23 L is stronger than that on 5 23 R . In Fig. 4(b) we show the 
bosonic u 2 decay branching ratio B(u 2 — > uih°) in the S 2 RL -5 3 RL plane. Its dependence on 
5 2 RL and 5 3 RL can be explained by that of \c{u 2 U\h°)\ 2 shown in Fig. 3(b). The branching 
ratio B{u 2 — > Hih ) can go up to 60% satisfying all relevant experimental and theoretical 
constraints listed in Appendix B. At the reference point of Table 1 it is about 46%. Note 
that at this point B(u 2 c%?) — 40% and B(u 2 tx\) — 10% whereas B(u 2 uiZ°) is 
very small (~ 1%). Note also that QFV effects (mainly due to the c R — t R mixing and the 
c R — ti mixing (~ Tu 32 )) together with the large top trilinear coupling Tu 33 (~ i R — II 
mixing parameter) can lead to a strongly enhanced branching ratio for the u 2 decay into 
h°. 

We have also studied the QFV decays u 3 — > Uih° and u 3 — >■ u 2 h°. The branching ratio 
of u 3 — > u\h° can go up to 30%. At the reference point of Table 1 is about 22%. u 3 has a 
large II component hence the u 3 — > uih° transitions are mainly due to II — > t R h°. On the 
other hand, in the decay u 3 — > u 2 h° the behaviour of the branching ratio is very different, 
because the II — c R transitions are more important. Its branching ratio at the reference 
point is about 10%. The branching ratios of u 3 — > uiZ° and u 3 — > u 2 Z° at the reference 
point are about 28% and 14%, respectively. 

In Scenario A the squark u 2 can also be produced in the decay of the gluino. We show 
in Fig. 5(a) and 5(b) the branching ratios of the decays g — > u 2 c + c.c. and g — > + c.c. 
as functions of S 2 RL and 5' 33 L . The branching ratio of the decay g — > u 2 i + c.c. is much 
smaller than that of g — > u 2 c + c.c. due to phase space. 

Furthermore, we have studied a scenario (Scenario B) where the gaugino mass param- 
eters Mi, M 2 and M 3 obey the GUT relation M 1 w 0.5 M 2 , M 3 /M 2 = g$/g%, where g 2 
and g 3 are the SU(2) and SU(3) gauge coupling constants, respectively. More definitely, 
we have changed only M 1 , M 2) and M 3 with respect to the values of Table 1, leaving all 
other parameters unchanged. We choose M x = 250 GeV, which leads to M 2 = 500 GeV 
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Figure 4: The branching ratio of the decay u 2 — » Wi/i as a function of $23 and <5. 
(a) and as a function of #23 and ^-f 1, (b) in Scenario A, with parameters as shown in 
Table 1. In both plots "X" indicates the reference point. 
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Figure 5: The branching ratios B(g — > u 2 c)+B(g — > u 2 c) (a) and B(g — > u 2 t)+B(g — >■ i^) 
(b) as functions of and 5^ L in Scenario A. The corresponding parameters are given 
in Table 1. "X" indicates the reference point in both plots. 



and M3 = 1500 GeV. The physical masses of the squarks are almost the same as in 
Table 2. In this scenario the gluino is relatively heavy, m g = 1626 GeV, therefore, it has 
a relatively small production cross section pp — > ggX. As we will see in the next section, 
gluino production is important, because the lighter squarks Ui >2 are also produced in the 
gluino decays g — > U\, 2 cj- Concerning the mass of the h°, the coupling |c(-UiW2^ )| 2 and 
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Figure 6: S^ L dependence of the branching ratios B(£t 2 — > U\h°) and B(tt 2 — y uiZ°) in 
Scenario C. "X" indicates the reference point defined with the parameters of Table 1, 
except for those shown in eq. (15). The vanishing of B(u 2 — > u\h°) at $J3 ~ —0.13 is 
due to kinematics, which also causes the peack of B(w 2 — > uiZ°). 

the branching ratio u 2 — y Uih°, they are very similar to those shown in Figs. 2, 3, and 4. 

In exploring various scenarios we have also studied one (Scenario C) where the QFV 
bosonic squark decays u 2 — y u\h° and u 2 U\Z° dominate over the fermionic u 2 decays. 
We have again changed some of the MSSM parameters with respect to Table 1, leaving 
all other parameters unchanged, 



In particular, the QFV trilinear coupling parameter 5 23 L {^ Tu3 2 ) is much larger than that 
in Scenario A. This new scenario satisfies all experimental and theoretical constraints listed 
in Appendix B. For the physical masses and branching ratios, we obtain the following 
results: m- g = 1134 GeV, m„ 1 = 651 GeV, m i2 = 800 GeV, = 1580 GeV, m^o 
= 125.3 GeV, m^= 398 GeV, B(u 2 -»■ ^h ) = 43%, B(u 2 u x Z°) = 34%, B(g -»■ 
u 2 c)+B(g -)• ii 2 c) = 8%, B(g -> u 2 t)+B{g ->• U 2 t) = 15%, B(6i -> £?c) = 96%, B(«i -)• 
X?i) = 4%. Note that both B(w 2 — )■ Wi/t ) and B(w 2 — )■ UiZ°) are very large, leading to 
the dominance of the QFV bosonic decays of u 2 . Note also that ui j2 are mixtures of cr 
and ti due to the sizable QFV trilinear coupling T[/ 32 (= (c R — £/,— iJ 2 coupling) ~ #23^)) 
which significantly enhances the QFV decay u 2 — y uih°. The large B(w 2 — y u\Z°) is 
mainly due to the sizable component in u± j2 in this scenario whereas in Scenario A the 
ti component is small. In Fig. 6 we show B(u 2 —y U\Z°) and B(£t 2 —y Uih°) as functions 
of 6%2 L . We see that the ratio B(u 2 — y U\Z°) /B(u 2 —> Uih°) is sensitive to the QFV as 
well as to the QFC SUSY parameters. 




= (650 GeV) 2 , M 2 33 = (1600 GeV) 2 , M 2 Q33 = 
0, 6^ R = 0, 5 2 t L = -0.17, ^ = -0.3. 



(780 GeV) 2 , 



(15) 
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4 Signatures 



In this section we will work out the signatures to be expected from the QFV decays of u 2 
into the lightest Higgs boson h° within the scenarios considered. The lighter squark states 
can be produced directly, pp — > UiUiX, pp — > u 2 u 2 X, or via gluino production, pp — > ggX, 
where at least one of the gluino decays into U\ or u 2 , g — > Ui,2 u ] ^1,2 c; Ui,2 t, with q being 
a light quark. In Table 3 we list the QFV processes leading to at least one Higgs boson h° 
in the final state. We assume that the t-quark can be identified, but the c-quark not. Of 
course, additional c-tagging would be very helpful. In the following t denotes top-quark or 
anti-top-quark. What concerns the final states from u 2 u 2 pah production, the final states 
with one t are explicit QFV signals whereas those with no t and 2t are apparent QFC 
signals. Concerning the final states from gluino pair production, those with one t and 3t 
are explicit QFV signals whereas those with no t and 4t are apparent QFC signals; those 
with 2t are explicit QFV (apparent QFC) signals in case they are "top-quark top-quark" 
("top-quark anti-top-quark"). Note here that the final states with "top-quark anti-top- 
quark" can stem from both QFV and QFC gluino decays, e.g. gg — > (cih°Xi) + (tch°Xi) 
and gg — > (tih°Xi) + ( c c^°X?)- Note also that the signal events with tt(or tt)jj, such as 
tt (or tt)jjh°IpTX , can practically not be produced in the QFC MSSM (nor in the SM). 
Because of the tiny SM background the detection of such di-top-quark signal events could 
be useful for discriminating between the QFC MSSM and QFV MSSM. 

At yfs = 14 TeV, for Scenario A, the production cross section for pp — >■ ggX is 
148 fb including one-loop SUSY-QCD corrections. For that we used Prospino2 [29] as 
this cross section is only very weakly dependent on the QFV parameters. The cross 
sections for pp — > u 2 u 2 X, pp — » gu\X, and pp — > gu 2 X are at tree-level 10 fb, 1 fb, and 

1.4 fb, respectively. For the calculation of these cross sections we have used FeynArts and 
FormCalc [30, 31]. All numbers for the signal cross sections given in this section include 
the charge conjugate final states. 

In Scenario A, at the reference point (see Table 1), the branching ratio for u 2 — > u\ 
is about 47% and that for u 2 — > uiZ° is about 1%. We also have B(ui — > cxi) = 36%, 
B(mi txi) = 64%, B(w 2 ->■ cxi) = 43%, B(u 2 txi) = 9%. Therefore the produced 
u 2 u 2 state goes with a probability of roughly 15% into the final state 2j + h° + tf T - In our 
scenario, the corresponding cross section for pp — > u 2 u 2 X — > 2j + h° + 1£t + X is almost 

1.5 fb. Note, however, that this signature can also occur in the QFC bosonic decays. 
On the other hand, the process pp — > u 2 u 2 X — > j + t + h° + J^t + X is QFV and the 
corresponding cross section is almost 2.8 fb. Even the cross section for pp — > u 2 u 2 X — > 
j + t + 2h° + $ T + X is about 1 fb. As the ratio B(u 2 UiZ )/ B(u 2 ->• Uxh ) is only about 
0.02 at the reference point, the signature u 2 u 2 — > 2j + Z° + h° + Jfi T has a probability of 
only 0.1%. 

The squarks u 2 are also produced in gluino decays. The branching ratios for the 
gluino decays at the reference point of Scenario A are B(^ — > u\c + c.c.) = 18%, B((? — >• 
Uii + c.c.)=55%, B(g — > u 2 c + c.c.)=24% and B(g — > u 2 t + c.c.)=3%. For the process 
gg — > u 2 u 2 jj — > 4j + h° + $t + X one gets, therefore, a probability of 0.8%, leading 
to a cross section for pp — > ggX — > u 2 u 2 jjX — >■ 4j + h° + tf T + X of about 1.2 fb. In 
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Table 3: Possible signatures expected from the QFV decays of u 2 into the lightest Higgs 
boson h°. t denotes top-quark or anti-top-quark; j denotes a c-quark jet. We also give 
the cross section to be expected in Scenario A, in case it exceeds 1 fb. 



processes 


final states containing h° 


pp ->■ u 2 u 2 X 


2j + h° + ^ T + X (1.4 fb) 
j + t + h° + $ T + X (2.8 fb) 
j + t + 2h° + # T + X (1 fb) 
2j + Z° + h° + $ T + X 
2j + 2h° + $ T + X 
2t + h° + # T + X 
2t + 2h° + # T + X 



processes 


final states containing h° 


pp ->■ ggX 


Aj + h° + $ T + X (2 fb) 
3j + * + /i° + # T + X (8 fb) 
2j + 2t + h° + $ T + X (4 fb) 
4j + Z° + h° + tf T + X 

Aj +2h° + # T + X 
3j + t + 2h° + $ T + X 
2j + 2t + 2h° + tf T + X 
3j + t + Z° + h° + tf T + X 
2j + 2t + Z° + h° + tf T + X 
At + h /Z (+h°/Z ) + # T + X 



the process pp —¥ ggX — > u 2 u 2 jjX the final state Aj + 2h° + Ifi T + X is possible with a 
probability of approximately 0.3%. The cross section for pp — > ggX — >■ 4j + 2/i° + ^ T + X 
is 0.5 fb. 

A further interesting process is pp — > ggX — > itiu 2 jjX — > Aj + h° + J^t + X, having a 
probability of 0.5%, giving a cross section of 0.8 fb. Therefore, one has a cross section of 
2 fb altogether for ->• ggX ->• 4j + /i° + $ T + X. The QFV final state 3j+t + h° + # T +X 
coming from pp — > ggX has a cross section of 8 fb. Correspondingly, the final state 
2j + 2t + h° + + X from pp — >■ ggX has a cross section of 4 fb and that of pp — )■ ggX — > 
3j + t + 2h° + lfi T + X is almost 0.9 fb. We give in Table 3 the cross sections corresponding 
to the final states as discussed. All other final states shown have cross sections smaller 
than 1 fb. 

Summing up the cross sections for all QFV final states with at least one h° one gets 
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Figure 7: Signal rates for the QFV bosonic squark decay signatures from gluino pair 
production (a) in Scenario A and (b) in Scenario C as functions of S^ 1 . The red solid 
(dashed) line corresponds to the pure QFV one t final state 3j + t + h° (Z°) + Jfix + X. 
The green solid (dashed) line corresponds to the QFV signal events with 2t coming from 
tt/tt/ti+2j + h° (Z°) + IpT + X. The blue solid line corresponds to the QFC signal events 
t + t + 2j + h° + JpT + X. The violet dashed line corresponds to the pure QFV one t 
final state 3j + 1 + h° + Z° + IpT + X. Note that the final states containing Z° are not so 
important in Scenario A and therefore not shown in (a). "X" indicates the corresponding 
scenario's reference point: for Scenario A defined with the parameters of Table 1, and for 
Scenario C defined with the parameters of Table 1, except for those shown in eq. (15). 



19 fb in the scenario considered. This means that one could expect about 1900 of such 
events assuming an integrated luminosity of 100 fb -1 . 

In Fig. 7(a) we show the signal cross sections for pp — > ggX — > 3j + t + h° + ]^T + X and 
pp ggX 2j + 2t + h° + Ipx+X in Scenario A as a function of S^ 1 - The red solid line 
corresponds to the pure QFV one * final state 3j + t + h° + $ T + X. The green solid line 
corresponds to the pure QFV 2t final state tt/tt + 2j + h° + Ifir + X plus the 2t final state 
signal coming from QFV gluino decays, tt + 2j + h° + J^t + X. Note that the contribution 
of the tt/tt final state events is exactly equal to the contribution of the QFV tt final state 
events due to the Majorana nature of the gluino. The blue solid line corresponds to the 
QFC signal events t + i + 2j + h° + tfx + X coming from the QFC gluino decays. The 
behaviour at 5%z L = 0.074 is due to the fact that the decay Wi — > t\\ is kinematically not 
possible whereas the decay — > tx± is still allowed (with 18% branching ratio). Note that 
the QFV cross sections do not vanish for S^ 1 = because the other QFV parameters 
5^ L and J^P* are not zero. — In the GUT inspired scenario (Scenario B) the gluino is 
much heavier and therefore the cross section a(pp ggX) is much smaller being 3.5 fb. 
The final state coming from pp — > U2U2X have the same cross section as in scenario A 
whereas the final states due to pp — > ggX are about a factor of 40 smaller. 

The third scenario (Scenario C) is characterized by a higher branching ratio B(tt2 - > 
U\Z°) = 34%. Therefore, one expects final states with Z° and h°. As the gluino mass is 
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very close to that of Scenario A, the cross section for pp — >• ggX is 148 fb. Consequently, 
the cross section a(pp — > ggX — > 3j + t + h° + Ipr + A) is 8.5 fb and a(pp — > ggX — > 
3j + t + Z° + IpT + A) is 6.8 fb. In Fig. 7(b) we show the cross sections analogous to 
those shown in Fig. 7(a), but for Scenario C. In addition, we also show the cross sections 
of the QFV final states containing Z°. The green dashed line corresponds to the pure 
QFV 2t final state tt/tt + 2j + Z° + Jf, T + X plus the 2t final state signal coming from 
QFV gluino decays, tt + 2j + Z° + Ip T + X. The contribution of the tt/tt final state events 
is again equal to the contribution of the QFV tt final state events. The violet dashed 
line corresponds to pure QFV one t final state 3j + t + h° + Z° + Ifir + A. Fig. 7(b) is 
symmetric in S^ L —> S^ 1 - 

We want to comment shortly on the background coming from SM processes and MSSM 
processes. In the SM the main background is the production of a Higgs boson H in 
association with top quarks, pp — > tiHX, where H is radiated off from top or antitop. The 
cross section at y^s = 14 TeV is about 400 fb. In these events, however, there is no missing 
energy, ftp. Further Higgs boson production processes are pp — > Z°Z°H; W + W~H. They 
will of course, constitute a background to the h° + jets + J^t- However, these processes 
do not contain a top in the final state. As discussed above, events with top (antitop) in 
the final states together with a h° are the most significant signatures for QFV. Single H 
production via top loops in gluon-gluon fusion as well as pp — > bbHX do not contain a 
top quark in the final state either. 

As we have shown, the final state with t + h° + jets has the largest cross section 
in the scenarios considered. Concerning the background within the general MSSM, the 
situation can be more complex. In the scenarios considered the charginos and neutralinos 
are relatively heavy, so that the decays of the lightest squarks into these play a minor 
role, except those into the lightest neutralino. If this is not the case the QFV signals will 
be less pronounced. Clearly MC studies including detector effects and appropriate cuts 
would be necessary to select the proposed QFV signals. However, this is beyond the scope 
of this paper. 

5 Summary 

In this paper we have studied the effects of QFV in the bosonic squark decays U2 — > 
Uih°/Z° and u 3 — > iti^hP/Z at the LHC. We have assumed mixing between the second 
and third squark generations, that is Cl,r — Il,r mixing. In our calculations we have 
taken into account all experimental constraints from B meson data on AMb s , B(6 — > s'y), 
B(B S fifi), limits on the gluino and squark masses, the latest data on the lightest Higgs 
boson mass, m h o ps 126 GeV, and the theoretical constraints on the trilinear couplings 
from the vacuum stability conditions. We have found that the branching ratio B(-£t2 — > 
uih°) can be larger than in the QFC case, and can go up to 60%, and that the decay 
v,2 — > u\h° gives access to the QFV trilinear couplings T[/ 32 and 71/23 • We have studied 
the signatures expected from the QFV decay u 2 — > U\h° at LHC with y^i = 14 TeV in 
three different scenarios. We have considered direct u 2 production pp — > U2U 2 X as well 
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as u 2 production in g decays via pp — > ggX. In two scenarios (A and C) we have taken 
rrig ~ 1100 GeV and m fi2 m 860 GeV and in the third scenario (B) m g 1600 GeV. 
The most pronounced QFV signature is 3j ; + t + h° + Ifir + -X", for example coming from 
pp — > ggX — > ui^iu2cX — > ui^tuihPcX — > cicch°]^TX , which can have a cross section 
up to 8 fb in Scenario A. A further interesting QFV signature coming from gluino pair 
production is 2j + 2t + h° + If T + X with a cross section in Scenario A of about 4 fb. 

In conclusion, our analysis suggests that for a complete determination of the MSSM 
parameters it will be necessary to study both the fermionic and the bosonic QFC and 
QFV decays of squarks. This can also have an influence on the squark and gluino searches 
at LHC. 
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A Up-squark decay into h° 



In the super-CKM basis, the Lagrangian including the coupling of up-type squarks to the 
lighter neutral Higgs boson, h°, is given by 



C = 



92 h° 



2m 



w 



/ \ cc 

u* L Uj L \m 2 w sin(a + (3)(1 - - tan 2 6 W ) 5^ + 2 — 

V o SI 



cos a: 



■ 3 m u i $ij 

sm (3 ' J 



+ u m u jR ( m VK sin (« + 0) \ tan2 Ow&ij + 2 ^^ m «.< <J «) 



+ 



_ _ / ,sma . cos a v 2 
u iR u jL (l* ^Tn m uAj + 



sm, 



sin/3 y/2 



a + h.c. 



(16) 



The terms going with m 2 ^ stem from the D-terms of the scalar potential and the 
expressions with quark masses m u ^ stem from Yukawa and F-terms. They are all flavour- 
universal. The trilinear couplings are explicit breaking terms that couple left-handed to 
right-handed squarks. Inserting the transformations to the physical fields, u iL = (R uji ) ik u k 
and itiR = (R uJ< )(i+3)kUk, eq. (16) can be written in terms of physical up-type squark fields 
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as C = Cu.u- h o u* Ui h° with the coupling 



92 



m 2 w sin(a + 0) \{l - \ tan 2 d w ){Ku) jk {1l\) kl 



2m w 

+ 1 tan 2 6 w (1Zu)j (fc+3) (Hu) (k+3) 



cos a 

+ 2 — ~R 
sinp 

sin a 



+ 



+ 



sin/3 
cos a V2 



{Ku)jk ml k {n\) ki + {nu) j{ k+z)rn\ k {n\) {k+z) 

j (k+3)VT7"u,k 



(72«)j(fc+3) (Tu)ik {K\)u + (^)«fc (^. fi )(i+3) 



sin/3 a/2 

The decay width for the process ui — > Ujh° is given by 



1 /<™L™L^°) 



16tt 2m^ 



As usual, k is defined by k 2 (x, y, = (x — y — z) 2 — 4yz. 



(17) 



(18) 



B Experimental and theoretical constraints 

Here we summarize the experimental and theoretical constraints taken into account in the 
present paper. The constraints on the MSSM parameters from the B-physics experiments 
and from the Higgs boson search at LHC are shown in Table 4. 

The particle discovered most recently at LHC [25, 32, 33] and consistent with the SM 
Higgs boson we identify as the MSSM Higgs boson h° which is indeed SM-like in the 
decoupling Higgs scenarios considered in our paper. For the mass of the Higgs boson h°, 
we take an average of the central values ((125.3(CMS) + 126.5(ATLAS))/2 ~ 126 GeV) 
and adding the theoretical uncertainty of ~ ±2 GeV [34] linearly to the experimental 
uncertainty of the CMS data at 1.5 a, we take 123 GeV < m^o < 129 GeV. Very 
recently LHCb has found a first evidence for B s — > at 3.5 a level and obtained 

B(B S — > /i + /U~) = (3.2 + 1.5/ — 1.2) x 1CT 9 , which is in very good agreement with the 
SM prediction [35]. They also set a double sided limit at 95% CL: 1.1 x 1(T 9 < B{B S ->■ 
yU+yU - ) < 6.4 x 1CT 9 . We have confirmed that this limit is satisfied in all plots presented in 
this article: for this branching ratio we have obtained numerical results very close to the 
SM prediction, which is mainly due to the fact that we take a very heavy Higgs mass m^o 
(= 1500 GeV) and that we assume almost no s — b mixing. We do not take the constraint 
from the recent Babar data on B(B — > D^tv) [36] because, as pointed out in [37], it 
is probable that the SM predictions are significantly underestimated and hence that the 
constraint on tan/3/m#+ would actually be much weaker than that claimed in [36]. 
In addition to these constraints we also impose the following experimental constraints: 
(i) The LHC limits on the squark and gluino masses (at 95% CL) [48, 49] 
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Table 4: Constraints on the MSSM parameters from the B-physics experiments relevant 
mainly for the mixing between the second and the third generations of squarks and from 
the limit on the h° mass. The fourth column shows constraints at 95% CL obtained by 
combining the experimental error quadratically with the theoretical uncertainty, except 
for B(B S w) and m h o. R^= B ff B V_I T + + ts" ~ [!-((™b+ tan /3)/m„+) 2 ] 2 where m H+ 
is the H + mass [38] . 



Observable 


Exp. data 


Theor. uncertainty 


Constr. (95%CL) 


AM Bs [ps- 1 ] 
10 4 xB(6^ s 7 ) 
10 6 xB(6->- s l+l~) 

(/ = e or /i) 
10 9 xB(B s -> 
10 4 xB(B+ -> t+i/) 

to/jo [GeV] 


17.725 ± 0.049 (68% CL) [39] 
3.37 ± 0.23 (68% CL) [41] 
1.60 ± 0.50 (68% CL) [43] 

< 4.2 (95% CL) [45] 
1.15 ± 0.23 (68% CL) [46] 

125.3 ±0.64 (68% CL)(CMS), 
126.5(ATLAS) [25] 


±3.3 (95% CL) [40] 
±0.23 (68% CL) [42] 
±0.11 (68% CL) [44] 

±0.25 (68% CL) [47] 
±2 [34] 


17.73 ± 3.30 
3.37 ±0.64 
1.60 ± 1.00 

< 4.2 

nSUSY _ 

B — >TV _ 

0.96 ±0.54 
123 < m h o < 129 



• In a simplified model: m g > 1050 GeV, m g > 1380 GeV (for m^o =0), no limit 
on rrig and rriq (for m^o > 400 GeV); 

• In the context of the CMSSM/mSUGRA: m g > 840 GeV, m- q > 1300 GeV. 

Here rriq represents the mass of mass-degenerate 1 st and 2 nd generation squarks. Note 
that in general, top-squark and bottom-squark mass limits are significantly weak and 
quite assumption dependent [48, 50]. Taking into account these limits, we assume a 
gluino mass larger than 1 TeV in our analysis. 

(ii) ATLAS and CMS set limits on m-± and m^o from negative searches for charginos 
and neutralinos in leptonic final states [51]. Respecting these limits we impose the 
following conditions: m~± > 400 GeV, m^o > 200 GeV. 

(iii) The constraint on (m^o , tan/3) from the MSSM Higgs boson searches at LHC 
[52]. 

(iv) The experimental limit on SUSY contributions on the electroweak p parameter 
[53]: Ap (SUSY) < 0.0012. 

Furthermore, we impose the following theoretical constraints from the vacuum stability 
conditions for the trilinear coupling matrices [54]: 

\T Uaa \ 2 < 3 Y* a (M 2 Qaa + M 2 Uaa + mj) , (19) 

\T Daa \ 2 < 3 Y 2 a (M 2 Qaa + M 2 Daa + mj) , (20) 

|7W| 2 < (M 2 Qaa + M 2 ^ + m 2 2 ) , (21) 

\T Da p\ 2 < Y 2 Dl (M 2 Qaa + M 2 D n + mj) , (22) 
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where a, j3 — 1, 2, 3, a j£ (3; 7 = Max(a:, /3) and m 2 = (m^ ± + m| sin 2 6*vk) sin 2 P — 
m 2 = (ra^± + m| sin 2 6*^ ) cos 2 f3 — \m 2 z . The Yukawa couplings of the up-type and down- 
type quarks are Y Ua = y/2m Ua /v 2 = ^ m "" s f n/3 («« = u,c,t) and Y Da = y/2m da /v 1 = 
^2 mwcosi3 (d a = d, s, b) , with m Ua and being the running quark masses at the weak 
scale and g being the SU(2) gauge coupling. All soft-SUSY-breaking parameters are also 
assumed to be given at the weak scale. As SM parameters we take my/ = 80.4 GeV, 
m z = 91.2 GeV and the on-shell top-quark mass m t = 173.3 GeV [55]. We have found 
that our results shown are fairly insensitive to the precise value of m t . 

We take into account all constraints listed in this section in all plots presented in this 
article. The constraints on the QFV parameters from B(b — > s'-f), AMb s , fn^o, and the 
vacuum stability conditions are especially important for this study. We use the public 
code SPheno v3.1 [22, 23] in the computation of the observables except for cross sections 
(i. e. physical masses, decay branching ratios, ^m Bs and A p (SUSY)) in this article. 

References 

[1] A. J. Buras et al, Phys. Lett. B 500 (2001) 161 [arXiv:hep-ph/0007085]. 

[2] G. D'Ambrosio, G. F. Giudice, G. Isidori and A. Strumia, Nucl. Phys. B 645 (2002) 
155 [arXiv:hep-ph/0207036]. 

[3] A. L. Kagan, G. Perez, T. Volansky and J. Zupan, Phys. Rev. D 80 (2009) 076002 
[arXiv:0903.1794, hep-ph]. 

[4] J. Beringer et al. [Particle Data Group Collaboration], Phys. Rev. D 86 (2012) 
010001. 

[5] G. Hiller and Y. Nir, JHEP 0803 (2008) 046 [arXiv:0802.0916 [hep-ph]]. 

[6] G. Hiller, J. S. Kim and H. Sedello, Phys. Rev. D 80 (2009) 115016 [arXiv:0910.2124 
[hep-ph]]. 

[7] M. Muhlleitner and E. Popenda, JHEP 1104 (2011) 095 [arXiv:1102.5712 [hep-ph]]. 

[8] G. Bozzi, B. Fuks, B. Herrmann and M. Klasen, Nucl. Phys. B 787 (2007) 1 
[arXiv: 0704. 1826 [hep-ph]]. 

[9] B. Fuks, B. Herrmann and M. Klasen, Nucl. Phys. B 810 (2009) 266 [arXiv:0808.1104 
[hep-ph]]. 

[10] A. Bartl et al, Phys. Lett. B 698 (2011) 380 [Erratum-ibid. B 700 (2011) 390] 
[arXiv: 1007.5483 [hep-ph]]. 

[11] M. Bruhnke, B. Herrmann and W. Porod, JHEP 1009 (2010) 006 [arXiv:1007.2100 
[hep-ph]]. 



18 



[12] T. Hurth and W. Porod, JHEP 0908 (2009) 087 [arXiv:0904.4574, hep-ph]. 

[13] A. Bartl et al, Phys. Lett. B 679 (2009) 260 [arXiv:0905.0132 [hep-ph]]. 

[14] A. Bartl et al, Phys. Rev. D 84 (2011) 115026 [arXiv: 1107.2775 [hep-ph]]. 

[15] B. Fuks, B. Herrmann and M. Klasen, Phys. Rev. D 86 (2012) 015002 
[arXiv: 11 12.4838 [hep-ph]]. 

[16] A. Bartl, W. Majerotto and W. Porod, Z. Phys. C 64 (1994) 499 [Erratum-ibid. C 
68 (1995) 518]. 

[17] A. Bartl et al, Phys. Lett. B 435 (1998) 118 [hep-ph/ 98042 65]. 

[18] A. Djouadi, J. L. Kneur and G. Moultaka, Nucl. Phys. B 569 (2000) 53 [hep- 
ph/9903218]. 

[19] B. C. Allanach et al., Comput. Phys. Commun. 180 (2009) 8 [arXiv:0801.0045 [hep- 
ph]]. 

[20] F. Gabbiani, E. Gabrielli, A. Masiero and L. Silvestrini, Nucl. Phys. B 477 (1996) 
321 [hep-ph/9604387]. 

[21] J. A. Aguilar-Saavedra et al, Eur. Phys. J. C 46 (2006) 43 [hep-ph/0511344]. 

[22] W. Porod, Comput. Phys. Commun. 153 (2003) 275 [hep-ph/0301101]. 

[23] W. Porod and F. Staub, Comput. Phys. Commun. 183 (2012) 2458 [arXiv: 1104. 1573 
[hep-ph]]. 

[24] F. Staub, T. Ohl, W. Porod and C. Speckner, Comput. Phys. Commun. 183 (2012) 
2165 [arXiv:l 109.5147 [hep-ph]]. 

[25] CMS data: J. Incandela, plenary talk at 36th International Conference on High 
Energy Physics, Melbourne, Australia, 4-11 July 2012; ATLAS data: R. Hawkings, 
plenary talk at 36th International Conference on High Energy Physics, Melbourne, 
Australia, 4-11 July 2012. 

[26] F. Brummer, S. Kraml and S. Kulkarni, JHEP 1208 (2012) 089 [arXiv: 1204.5977 
[hep-ph]]. 

[27] A. Arbey, M. Battaglia, A. Djouadi and F. Mahmoudi, JHEP 1209 (2012) 107 
[arXiv: 1207. 1348 [hep-ph]]. 

[28] R. Benbrik et al., arXiv: 1207. 1096 [hep-ph]. 

[29] W. Beenakker, R. Hopker, M. Spira, P. M. Zerwas, JHEP 08 (2010), 098 
[arXiv: 1006.4771 [hep-ph]]. 



19 



[30] T. Hahn, Comput. Phys. Commun. 140 (2001), 418 [arXiv:hep-ph/00 12260]. 

[31] T. Hahn and C. Schappacher, Comput. Phys. Commun. 143 (2002), 54 [arXiv:hep- 
ph/0105349]. 

[32] G. Aad et al. [ATLAS Collaboration], Phys. Lett. B 710 (2012) 49 [arXiv:1202.1408 
[hep-ex]]. 

[33] S. Chatrchyan et al. [CMS Collaboration], Phys. Lett. B 710 (2012) 26 
[arXiv: 1202. 1488 [hep-ex]]. 

[34] S. Heinemeyer, O. Stal and G. Weiglein, Phys. Lett. B 710 (2012) 201 
[arXiv: 11 12.3026 [hep-ph]]. 

[35] J. Albrecht, plenary talk at Hadron Collider Physics Symposium 2012 (HCP2012), 
Nov 12-16, 2012, Kyoto, Japan. 

[36] J. P. Lees et al. [BaBar Collaboration], Phys. Rev. Lett. 109 (2012) 101802 
[arXiv: 1205.5442 [hep-ex]]; G. D. Nardo, parallel talk at 36th International Con- 
ference on High Energy Physics, Melbourne, Australia, 4-11 July 2012. 

[37] M. Danilov, parallel talk at 36th International Conference on High Energy Physics, 
Melbourne, Australia, 4-11 July 2012. 

[38] W. -S. Hou, Phys. Rev. D 48 (1993) 2342. 

[39] S. Vecchi, parallel talk at 36th International Conference on High Energy Physics, 
Melbourne, Australia, 4-11 July 2012. 

[40] M. S. Carena et al, Phys. Rev. D 74 (2006) 015009 [hep-ph/0603106]; see also P. Ball 
and R. Fleischer, Eur. Phys. J. C 48 (2006) 413 [hep-ph/0604249]. 

[41] S. Stone, plenary talk at 36th International Conference on High Energy Physics, 
Melbourne, Australia, 4-11 July 2012. 

[42] M. Misiak et al, Phys. Rev. Lett. 98 (2007) 022002 [hep-ph/0609232]. 

[43] M. Iwasaki et al. [Belle Collaboration], Phys. Rev. D 72 (2005) 092005 [hep- 
ex/0503044]; B. Aubert et al. [BABAR Collaboration], Phys. Rev. Lett. 93 (2004) 
081802 [hep-ex/0404006]. 

[44] T. Hubcr, T. Hurth and E. Lunghi, Nucl. Phys. B 802 (2008) 40 [arXiv:0712.3009 
[hep-ph]]. 

[45] M. Perrin-Terrin, parallel talk at 36th International Conference on High Energy 
Physics, Melbourne, Australia, 4-11 July 2012. 

[46] R. Barbieri, summary talk at 36th International Conference on High Energy Physics, 
Melbourne, Australia, 4-11 July 2012. 



20 



[47] K. Trabelsi, plenary talk at ICHEP2010, PoS (ICHEP 2010) 566. 

[48] A. Parker, plenary talk at 36th International Conference on High Energy Physics, 
Melbourne, Australia, 4-11 July 2012; M. Backes, parallel talk at the same conference. 

[49] G. Aad et al. [ATLAS Collaboration], arXiv: 1208.0949 [hep-ex] (submitted to Phys. 
Rev. D); G. Aad et al. [ATLAS Collaboration], arXiv: 1208.4688 [hep-ex] (sub- 
mitted to Phys. Rev. D); G. Aad et al. [ATLAS Collaboration], ATLAS-CONF- 
2012-103; G. Aad et al. [ATLAS Collaboration], arXiv: 1206. 1760 [hep-ex] (sub- 
mitted to JHEP); S. Chatrchyan et al. [CMS Collaboration], arXiv: 1206.3949 
[hep-ex] (submitted to Phys. Lett. B); S. Chatrchyan et al. [CMS Collaboration], 
http://cdsweb.cern.ch/record/1460434; S. Chatrchyan et al. [CMS Collaboration], 
http://cdsweb.cern.ch/record/1460433; S. Chatrchyan et al. [CMS Collaboration], 
arXiv: 1208.4859 [hep-ex] (submitted to Phys. Rev. D); S. Chatrchyan et al. [CMS 
Collaboration], arXiv:1205.6615 [hep-ex] (submitted to Physical Review Letters); S. 
Chatrchyan et al. [CMS Collaboration], arXiv: 1204.3774 [hep-ex] (submitted to 
Physics Letters B). 

[50] A. Cakir [CMS Collaboration], arXiv: 121 1.6289 [hep-ex]; G. Aad et al. [ATLAS 
Collaboration], arXiv: 1209.2102 [hep-ex] (submitted to Phys. Lett. B); G. Aad et al. 
[ATLAS Collaboration], arXiv: 1208.4305 [hep-ex] (submitted to Eur. Phys. J. C); G. 
Aad et al. [ATLAS Collaboration], arXiv: 1208.2590 [hep-ex]; G. Aad et al. [ATLAS 
Collaboration], arXiv: 1208. 1447 [hep-ex] (submitted to Physical Review Letters); G. 
Aad et al. [ATLAS Collaboration], arXiv: 1207.4686 [hep-ex] (submitted to Eur. Phys. 
J. C); G. Aad et al. [ATLAS Collaboration], arXiv: 1203.6193 [hep-ex] (submitted to 
Physical Review D); G. Aad et al. [ATLAS Collaboration], arXiv:1203.5763 [hep-ex] 
(submitted to Physical Review Letters); G. Aad et al. [ATLAS Collaboration], Phys. 
Rev. Lett. 108 (2012) 181802 [arXiv: 1112.3832 [hep-ex]]. 

[51] G. Aad et al. [ATLAS Collaboration], arXiv: 1208.3144 [hep-ex]; S. Chatrchyan et al. 
[CMS Collaboration], [arXiv:1209.6620 [hep-ex]]. 

[52] S. Chatrchyan et al. [CMS Collaboration], Phys. Lett. B 713 (2012) 68 
[arXiv: 1202.4083 [hep-ex]]. 

[53] G. Altarelli, R. Barbieri and F. Caravaglios, Int. J. Mod. Phys. A 13 (1998) 1031 
[hep-ph/9712368]. 

[54] J. A. Casas and S. Dimopoulos, Phys. Lett. B 387 (1996) 107 [hep-ph/9606237]. 
[55] E. Shabalina, plenary talk at ICHEP2010, PoS (ICHEP 2010) 561. 



21 



